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I. INTRODUCTION 



It was Ehrenfest (Enzykl. Math. Wissensch, IV, 2(11) fasc.6, p82, note23, 1912) who first 
asked whether a function exists which, like entropy in the equilibrium state of an isolated 
system, achieves its extreme value in a stationary non-equilibrium state. 

There are various results and formulations of irreversible thermodynamics based on the 
extremum of entropy production. They are related either to the minimum or maximum 
entropy production. The well known example based on minimum entropy production is 
Prigogine theorem [l|. 

There are several different results that concern maximum entropy production (MEP). 
Ziegler applied the principle of maximum entropy production in thermomechanics [2]. It 
states that the rate of entropy production under prescribed forces should be maximum. Pal- 
tridge {3] has formulated empirically MEP principle in order to describe the Earth climate. 
Apart from details of his model (like albedo or cloudiness) the main features of his model are 
steady incoming sun radiation, outgoing radiation proportional to the fourth power of the 
local annual average temperature and horizontal heat flux from equator to poles. There are 
no boundary conditions on atmosphere. Paltridge has proposed that steady state is the state 
of maximum entropy production due to the latitudinal heat transport. His predictions of the 
annually average distribution of the temperature of atmosphere ( and some other climatic 
parameters) fit very well measured values. Kohler {4] has started from the Boltzmann's 
transport equation and has shown that the solution for the velocity distribution function of 
rarefied gas, in stationary state close to the equilibrium, is the state of extremum entropy 
production. He has argued that extremum type depends on the choice of the constraint. 
We show in this paper that constraint that leads to the stationary state as the state of the 
minimum entropy production is not consistent with his starting assumption. We argue here 
that rarefied gas in the stationary state close to equilibrium is in the MEP state. 

Recently Jaynes' principle of maximum information entropy (MaxEnt) has been exploited 
to derive the MEP principle. Dewar has introduced paths as possible trajectories of 

the system in phase space. Then using the MaxEnt procedure he has shown , under some 
assumptions, that the most probable development of the system is accompanied with MEP. 
Niven [?| has applied MaxEnt using the values of fluxes as variables and he found that the 
stationary state of the system is the MEP state. 
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At first sight minimum and maximum entropy production results seem to contradict 
each other. But if one considers starting assumptions one finds that these assumptions are 
different. Thus these results are independent of one another j^l, 

In this paper we focus on Onsager's principle of least dissipation of energy and show that 
this principle is equivalent to MEP principle. Secondly, as we have already noted, we show 
that Kohler interpretation of stationary state of rarefied gas close to the equilibrium should 
be restricted only on the constraint that it is consistent with the interpretation of this state 
as the MEP state. 

This is the first of three papers in a series and are referred to as papers I, II and III. In 
paper II we discuss whether stationary or relaxation processes are suitable for the formulation 
of principles. We argue in favour of relaxation processes. In paper III we apply MaxEnt 
formalism to relaxation processes and we derive the MEP principle as its corollary. 



II. THE PRINCIPLE OF THE LEAST DISSIPATION OF ENERGY AND LINEAR 
NONEQUILIBRIUM THERMODYNAMICS 



Onsager's famous papers [lOi , examine linear nonequilibrium thermodynamics. The 
first sections of these papers examine reciprocal relations. The following sections are de- 
voted to the formulation of linear nonequilibrium thermodynamics. Onsager has employed 
;he principle of the least dissipation of energy. In the second part of his second paper 



11] , Onsager has applied the principle of the least dissipation energy to the general linear 
nonequilibrium process. 

A linear relationship between n independent thermodynamic forces {X{\ and their con- 
jugated fluxes {ji}, 

3 

exists for the system close to equilibrium. The density of the entropy production is equal to 
the product of the thermodynamic forces and conjugated fluxes 12]. The total entropy 
production is 



/ /(, 
E X ^dV. (2) 



a 
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In his second famous paper Onsager has introduced the dissipation function 



i r n 

~ 2 Yl R ijJdjdV, (3) 

J i,j=l 

and formulated the principle of the least dissipation of energy, 

(r — $ = maximum. (4) 

The variational procedure 

5 [a - $] = 0, (5) 

gives the system of equations ([1]). In other words, the linear irreversible thermodynamics 
can be inferred from the principle of the least dissipation of energy. 

III. THE MEP PRINCIPLE AND LINEAR NONEQUILIBRIUM THERMODY- 
NAMICS 

An alternative approach to linear nonequilibrium thermodynamics is based on the phe- 
nomenological fact that nonequilibrium processes are characterized by fluxes. Therefore, 
physical quantities relevant for a description of the time development of the system must be 
functions of fluxes. The standard approach to the nonequilibrium thermodynamics is based 
on the laws of conservation of mechanical physical quantities mass, momentum and energy. 
The next step is to add heat as the additional mechanism of the exchange of energy between 
systems. Then assuming local equilibrium one comes to the density of entropy production 
written as the sum of products of heat flux and thermal thermodynamic force and viscous 
pressure tensor and strain tensor divided by temperature. This result serves as the basis 
for the canonical form of the second postulate of irreversible thermodynamics that reads: 



Entropy production can be a 



ways written as the sum of products of thermodynamic forces 



and corresponding fluxes 12|, [13|]. Thermodynamic forces and fluxes are not independent 
quantities. In principle fluxes determine forces and vice versa. Having in mind this fact and 
the second postulate of irreversible thermodynamics we can say that the dynamical state of 
the system is described only by fluxes, or equivalently by thermodynamic forces. 

Entropy production is a basic, characteristic quantity of a nonequilibrium state. If the 
system is close to the equilibrium state we can make the Taylor expansion of the density of 
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entropy production da/dV in fluxes up to the second order 

da 
dV 



% = A + H B ^ + H R ddr (6) 



■th 



Here ji is the mean value of the i flux. The first term on the right-hand side of this equation 
vanishes since there is no entropy production in the equilibrium state. Coefficients Bi and 
Rij are the property of the system in equilibrium state. Entropy production does not depend 
on the direction of the flux flow, that is to say, it must be invariant under the replacement 
of {ji} {—ji}- This means that coefficients that multiply the odd power of fluxes vanish, 
that is, Bi = 0. In the lowest order, entropy production is the bilinear function of fluxes, 



a 



Y\ / RijJiJjdr > 0. (7) 
" " Jv 



Comparison of Eqs. and ([2]) shows that the dissipation function is in fact entropy 
production written in the space of fluxes. When a system is close to equilibrium with locally 
well defined intensive thermodynamic quantities the entropy production is given by Eq. ([2]), 
fl, El- 

It is pointed in the introductory part of this Section that canonical form of entropy 
production (|2]) comes from the law of conservation of energy. Then equation, 



/ XiJi ~ Ritfrii I dr = °- 



(8) 



is the constraint imposed the entropy production due to the law of conservation of energy. 
We seek the maximum entropy production (CO) taking into account the constraint (JED- 



The standard procedure 



14| is to find the maximum of the functional 



l i3 



XiJi ~ Hi J JJ > 



dr. (9) 



A standard variational calculus of extremum values combined with the constraint ([8]) gives 
A = 2 and Fi becomes 



F = -J2 J Rijjijjdr + 2 I J2 X didr. 

i,j V V i 



This is an equation of a quadratic surface turned upside down. The corresponding ex- 
tremum of F is maximum. The system of equations that determines the point of maximum 
, dF/dji = 0, is just a system of linearly coupled thermodynamic forces and fluxes ([!]). 



From Eq. (J7J) follows Rij = Rji i.e. Onsager's reciprocal relations. In short, linear 
nonequilibrium thermodynamics and Onsager's reciprocity relations follow from the MEP 
principle. 

A problem analogous to this is the problem of biochemical cycle kinetics close to the 
equilibrium state. Starting from the assumption that the energy conservation law is valid 
for a whole network of biochemical reactions we have showed that fluxes are distributed in 
such a way to produce maximum entropy [15j . 



In references [16|, [lTj we have considered a linear planar electric network held at a constant 
temperature. There is only a temperature gradient between network and surroundings. 
Physically this can be achieved using thermally high and poor conducting material for a 
network and surroundings, respectively. There is no coupling between electric and heat 
currents. Only electric currents are coupled via electromotive forces. Assuming that the 
stationary state is the state of maximum possible generated heat the Kirchhoff loop law is 
derived The principle of the maximum heat dissipation is closely related to the 

MEP principles. 

Using mesh currents the law of charge conservation has been taken implicitly. The energy 
conservation law is used explicitly as the constraint. If one only takes charge conservation 
aw one comes to the conclusion that stationary state is state of minimum generated heat 



18l | . We shell dwell more on this question in the paper II submitted to this special issue. 

We note that two of us (P.Z. and D.J.) have considered the heat flow in the anisotropic 
crystal. It is shown in this special example that the principle of the least dissipation of 



energy applied by Onsager in reference 10] is equivalent to the MEP principle 19] 



IV. THE LINEARIZED BOLTZMANN EQUATION AND THE EXTREMUM OF 
ENTROPY PRODUCTION 

The important influence of constraints at extrema of entropy production can be seen in 
the case of linearized Boltzmann equation. 



Here we follow the elegant approach given in reference js]. The Boltzmann equation is 
valid for rarefied gas where collisions are very well defined events. One seeks the one-particle 
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distribution function /(r, v, t) satisfying the equation, 

df df df Fdf 

777 = + V 7T + — 77~ = L 10 
at at or m ov 

Here the two terms on the left-hand side of equation describe the change in the number 
of molecules in a given element of the phase space due to the collisionless motion of the 
molecule in the outer field F. The right-hand side of the equation describes the net change 
in the number of molecules in a given element of the velocity space due to molecule collisions. 
Assuming the instantaneous change of molecule velocities in the collisions and ta king into 



account the conservation laws in collisions the integral can be written in the form 20] 




(ff-ff)qda(q,e)dV. (11) 



Here / = f(r,v,t) ; /' = /(r, v',£) and / = /(r, v, t) ;/' = /(r, v', £) are the distribution 
functions of particles before and after collision, respectively. In addition, q — v' — v is the 
magnitude of relative velocity of particles before collisions, e = (v — v)/\v — v\ is a unit 
vector parallel to the velocity change of one particle in collision and da(q, e) is the differential 
cross section. It is known from collision theory that differential cross section depends on 



interaction potential, relative velocity q and unit vector e, i.e. on the scattering angle 20]. 

In the case of local equilibrium, intensive variables (temperature, concentration) are well 
defined functions of the space. Then one assumes the approximate solution of Eq. ({TO]) Q], 



19 



9|. 

/ = /o(l + *(v)), (12) 

where 

/„ = ra(— — ) 3/2 exp( — 13 

JU K 27rkT J FK 2kT J v ; 

is the Maxwell-Boltzmann distribution. 

Due to the assumption of local equilibrium the perturbed function must not contribute 
to the prescribed intensive parameters like temperature, mean velocity and density. Then 

J f ^(v)dv = 0, J f o ^(v)Vdv=0, y"/ ^(v)^dv = 0. (14) 

Now the collision integral up to the first order of the perturbed term becomes 

J (/) =111 /o/o(* + *'-*- *')qd<r(q, e)dV = 6*. (15) 
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The designation of \1/ functions is the same as the distribution functions in Eq. (TTTj) . 
Operator O 0], jsj is linear and has the following properties 

0(aA + pB) = aO A + 0OB, (16) 



AOBdw = / BOAdv, (17) 



AOAdv > 0. (18) 
If we designate the left-hand side of Boltzmann equation (fiU]) with —Z it becomes 

Z=-OV. (19) 

Multiplying the linearized Boltzmann equation with —k In / and integrating over velocity 
space we get 

ds r 

+ v-j fl = k / ^O^dv. (20) 
<9t " J 

Here s = —k J f In /dv is the density of entropy and j s = —k J vfhifdv is the flux of 
entropy. Entropy looks like a fluid. The left-hand side of Eq. ( 1201) is the total time change 
of the density of entropy, i.e. it is the density of the entropy production. Then, according to 
the aforementioned equation collisions between molecules described by function \l/ are the 
sources of entropy production. There are a lot of functions obeying conditions (fl4j) . We 
choose those distribution functions that obey the additional condition 



YZdv= - J YOYdv. (21) 

Function Y is not the solution of the Boltzmann equation, i.e. 

Z ^ -OY. (22) 

Y can be interpreted as a trial distribution function. The corresponding entropy production 
due to this distribution is 

^ + V-(j s ) Y = k[Y,dY}. (23) 

If the temperature gradient in x direction is the only thermodynamic force, the left-hand 
side of the equation ( 12~TT) can be written as [4|] 



Here (q x )Y is the heat flux associated with the distribution function Y. This equation gives 
a physical meaning to condition (12 ip. Varied distribution functions are selected in such 
a way that the entropy production due to the molecular collisions is equal to the entropy 
production due to the heat conduction. In other words entropy produced by molecular 
collisions is equal to the product of the heat flux and conjugated thermodynamic force. 
This condition is in accordance with the starting assumption of a stationary process. 

The multiplication of the Boltzmann equation with the perturbative distribution function 
\l/ and integration over molecular velocity space v gives 

J ^!Zdw = - J md^dv. (25) 

Using conditions (I16H18P we get, 

y"(* - Y)6{^ - Y)dw > 0, (26) 



*()*r/v / YOYdv- I i'OYdv- I YO\V,lv= / *0*r/v+ / YOYdv-2 / l Y)Wv>0, 

(27) 

Otf dv+ / YOYdv+2 I YZdv = I ^O^dv+ I YOYdv-2 I YOYdv = I ^O^dv- I YOYdv > 0, 



(28) 

^O^rfv > J YOYdv. (29) 

The left-hand side is the entropy production of the rarefied gas. It comes from Eq. (129]) 
that solutions of the Boltzmann equation are in accordance with MEP principle. 



There is another approach to this problem which has been proposed by Kohler {4]. Due 
to the fact that Eq. ( I24p is proportional to the heat flux he had fixed the heat flux and 
varied the distribution function under constraint 

q x = constant. (30) 

We note that due to the already fixed temperature field the entropy production is also 
fixed. Kohler did not impose stationary conditions on the trial distribution function. In his 
approach the entropy generated by molecular collisions need not be equal to the entropy 
produced by fluxes and gradients. If one does not take care of this unconsistency the 
conclusion can be reached, by using the calculus of variation, that solutions of the Boltzmann 
equation are those that generate minimum entropy production 



V. CONCLUSION 



The Onsager principle of the least dissipation of energy is valid for processes close to the 
equilibrium state. In this paper the equivalence between the Onsager principle of the least 
dissipation of energy and the MEP principle is established. Starting from the fact that fluxes 
are the main phenomenological characteristic of irreversible processes we have expanded the 
density of entropy production as the function of fluxes up to the second order. We have 
found that the dissipation function introduced by Onsager is the entropy production in the 
space of the fluxes. Invoking the first law of thermodynamics the equivalency between the 
Onsager principle of the least dissipation energy and the MEP principle is established. 

It follows from the Boltzmann equation that entropy production is quantitatively related 
to the collision integral. Collisions between molecules in nonequilibrium state produce en- 
tropy. In this paper the solutions of the linearized Boltzmann equations are considered. It 
is found that these solutions correspond to the extremum of entropy production. The na- 
ture of the extremum depends on the constraints. Assuming that the entropy produced by 
molecular collisions is equal to the entropy productions due to the heat conductions or/and 
viscosity we find that the solutions of the Boltzmann equation satisfy the MEP principle. 

The principle of minimum entropy production is valid if one grants fixed fluxes. However, 
the starting assumption of the fixed thermodynamic forces and the additional assumption 
of fixed fluxes leave no room for variation. Although both approaches are equivalent from 
the mathematical point of view one has to notice that the assumption of fixed fluxes is not 
consistent with the starting assumption of the stationary process. Namely, the additional 
request that entropy produced by molecular collisions need not be equal to the entropy 
production described by thermodynamic forces contradicts the starting assumption of the 
stationary process. 

In short, revisiting the linear nonequilibrium thermodynamics and linearized Boltzmann 
equation shows that both approaches are in accordance with the MEP principle. Generally 
we can conclude that the MEP principle is valid for processes close to the equilibrium state. 
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